h corrections in semiclassical formulas for smooth chaotic dynamics 
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The vahdity of semiclassical expansions in the power of h for the quantum Green's function have 
been extensively tested for billiards systems, but in the case of chaotic dynamics with smooth 
potential, even if formula are existing, a quantitative comparison is still missing. In this paper, 
extending the theory developed by Gaspard et al., Adv. Chem. Phys. XC 105 (1995), based on the 
classical Green's functions, we present an efficient method allowing the calculation of h corrections 
for the propagator, the quantum Green's function, and their traces. Especially, we show that the 
previously published expressions for h corrections to the traces are incomplete. 
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I. INTRODUCTION 

Gutzwiller's work has now become a milestone in the 
understanding of the properties of a quantum system 
whose classical counterpart depicts chaotic dynamics [y . 
Starting from Fcynman's path formulation of quantum 
mechanics, he has been able to complete the early stud- 
ies of Van Vleck |2], deriving expressions for the semi- 
classical propagator, and from this, for the quantum 
level density: the well-known Gutzwiller trace formula. 
The later is an asymptotic series in h and can be sep- 
arated into two parts; the leading order corresponds to 
the Thomas-Fermi (or extended Thomas-Fermi when in- 
cluding h corrections) average density of states B; the 
other part corresponds to the oscillations around the pre- 
ceding term and involves contributions from all periodic 
orbits of the system. This formula has been widely used 
to obtain approximate values for the quantum energy 
eigenvalues of classically chaotic systems: the hydrogen 
atom in magnetic field [@, |5_^ the helium atom g, 0, B|, 

resonant tmmcl diode [9|, 
Since then, the Gutzwiller 
trace formula has also been generalized to take into ac- 
count contributions of other kinds: diffractive effects p4[ , 
continuous families of periodic orbits p^, p^ M, ghost 
orbits, etc. 

At the same time, because the trace formula as de- 
rived by Gutzwiller only contained the leading term of 
the asymptotic expansion of the quantum level density, 
the systematic expansion of the semiclassical propaga- 
tor in powers of h has been the purpose of several stud- 



anisotropic Kepler problem 
biUiards ficl, fll |l2|, |ll, etc 



mm 



|17| . However, these corrections to the trace 
formula have only been tested for billiards, for which both 
classical and quantum properties are easier to calculate. 
In the present paper, we will show how, for quantum 
systems whose Hamiltonian separates into kinetic and 
smooth potential energies, h corrections can be computed 



with great accuracy, extending the method described in 
Refs. pa, O, based on classical Green's functions. Espe- 
cially, we will show that the previous derivation ||l^, |l^ 
of the correction to Gutzwiller trace formula is partially 
wrong. 

From a numerical point of view, all quantities involved 
in the calculation of the h corrections for a given classical 
path can be obtained as solutions of sets of first order dif- 
ferential equations to be integrated along this path using 
standard time integrators like the Runge-Kutta method. 
The number of equations in these sets can be quite large 
and can be probably reduced with a deeper analysis of 
their structures, in the same way that the amplitude in 
the Gutzwiller trace formula for a two-dimensional (2D) 
system can be obtained by integrating only a (2 x 2) 
matrix and not the whole monodromy matrix [ |l8| . How- 
ever, it would give rise to more complicated expressions 
and probably to additional difficulties in the numerical 
implementation, whereas the expressions given in the pa- 
per can be put in the computer as they stand. Also, the 
amount of CPU time and the memory needed by the 
codes are small enough, so that, on a first stage, the re- 
duction of the number of equations can be skipped. 

The paper is divided as follows: in Sec. II, after a brief 
description of the derivation of h corrections for the prop- 
agator _?i'(q, qo, r), we explain how the classical Green's 
functions and h corrections can be efficiently computed. 
In Sec. pU we explain how to obtain expression of the 
additional terms, arising from the time to energy domain 
transformation, in h corrections for the quantum Green's 
function G(q, qo, T). In the case of the trace of the prop- 
agator, a detailed derivation of the h corrections is car- 
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ried out in Sec. ^, leading to the proper formulas, along 
with the way they can be computed. The time to energy 
transformation is explained in Sec. |y|, leading to the h 
correction expression in the case of the quantum Green's 
function. Finally, section VI shows how to apply theoret- 
ical expressions obtained in the four preceding sections 
in the case of the 2D hydrogen in magnetic field and 
emphasizes the excellent agreement with numerical coef- 
ficients extracted from exact quantum calculation, using 



harmonic inversion |19, 2^, 21 



II. THE PROPAGATOR A'(q,qo,T) 

A. Feynman path integral 

The starting point is the Feynman path integral, whose 
discrete version, for a time independent Hamiltonian 
which separates into kinetic and potential energies, H = 
p^/2 + y(q), reads as follows O: 



following equation jTs) : 

V-g{t,t') = t6{t-t'), 



K{(i, qo, T) == / dqi dq2 , • • • , dqN-iCiTrihAty 



Nf/2 



X exp 



n=0 



At 



q„ At + 0{At) 



(1) 



where At — T/N, qjv = q and i(q, q) is the classical 
Lagrangian. 

For small values of h (i.e., the semiclassical limit), using 
the stationary phase approximation, all preceding inte- 
grals are expanded around the stationary solutions, that 
is the classical orbits q^'(t) going from qo to q during time 
T, each of them thus giving a contribution Ki{q,c{Q,T) 
to the propagator, whose final expression reads formally 
as follows [ p^ : 

i^,(q,qo,T) = Af)(q,qo,T){l + iSCi(q,qo,T) + 0(ft2)}, 

(2) 

where Ki (q, qo,T) is the dominant semiclassical con- 
tribution to the propagator K{q,qo,T): 



40) 



i^r(q>qo,T) 



{2mh)fl^ 



det 



9qi9qo 



Wf{ci,cio,T) 



1/2 



X exp 



^<(q,qo,r)-*|^Z 



(3) 



where W^f'(q, qo,T) is the classical action and vi is the 
Morse index of the orbit. The Ci(q, qo, T) expression is 
given by [0|: 



dtvl'^i{t)g,j{t,t)gki{t,t) 



TpT 



+ dtdt'v}^{t)v,lZ{t') 



r{3) 



1 
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X [3g^j{t,t)gki{t,t')g„,^{t',t') 

+2guit,t')g,^{t,t')gkn{t,t')], (4) 

where the V^^' (t) are higher-order derivatives of the po- 
tential V, evaluated at qf{t). 

The classical Green's function g{t,t'), associated with 
the classical orbit, is an (/ x /) matrix solution of the 



(5) 



where V is the Jacobi-Hill operator, controlling the linear 
stability around the classical orbit in the configuration 
space 



13 



V ■ 



dt2 



dqdq 



[<f\t)] 



(6) 



Furthermore, the fact that both initial and final point are 
fixed in the propagator K{q,qo,T) imposes the follow- 
ing boundary conditions on the classical Green's func- 
tion ill: 



g{o,t') = g{T,t') = o vt'e[o,r] 



B. Classical Green's function 



(7) 



If qi(T) is a conjugate point of qo, then the determi- 
nant det (— 9qq W[^) in formula (0) is formally infinite, 
but this happens only for restricted values of T, so that, 
in this section, we will focus on the general case, for which 
q;(r) and qo are not conjugate points. 

Apart from t — t' , g{t,t') obeying the homogeneous 
Jacobi-Hill equation P • ^ = 0, so that, introducing the 
notations 



g^{t,t')=g{t,t') iorO<t<t', 

g+{t,t')^g{t,t') ioit'<t<T, 



one immediately obtains 



fG±{t,t') 



M{t) 



A±{.t') 
B±{t') 



(8) 



(9) 



where M{t) is the (2/x2/) monodromy matrix, depicting 
the linear stability around the classical orbit in the phase 
space. A± and B± are four (/ x /) matrices, whose 
values are determined from the boundary conditions at 
time t = t': 



g+{t',t')-g^{t',t') = o, 



(10) 



and at times t — and t ^ T: 



g-{o,t')^o, 
g+(r,t')-o. 



(11) 



For a Hamiltonian which separates between kinetic and 
potential energy H — p^/2 -I- V{q), M{t) has the follow- 
ing simple structure: 



Mit) = 



■hit) Ji{t) 

Mt) Mt) 



(12) 



which leads us to the following explicit expressions for 
the four matrices A± and B±: 




0, 



jJ{t')-J-\T)MT)J^{t'), 
-Jr\T)MT)J^{t'), 



(13) 



provided that Ji^{T) is invertible. Ji(T) being the up- 
per right (/ X /) submatrix of the matrix M, gives the 
linear displacement of the final position for a change in 
the initial momentum (the initial position being fixed 
to qo), i.e., Sq{T) = Ji{T)5pQ. Thus, Ji{T) is the in- 
verse matrix of {—d^^ W^;'^') which has been supposed to 
be invertible (q(T) and qo are not conjugate points). Fi- 
nally the full expression for the classical Green's function 
reads: 



G{t,t') 



{ Mt) [jJ{t')-J^\T)UT)jJ{t')] 
for < t < i', 

[J2{t)~h{t)J-\T)MT)\ Jj{t') 
for t' <t<T. 

(14) 



Using the symplectic structure of M{T), one can show 
that 



g{t',t) = g^{tX} 



(15) 



as expected because the operator V and the boundary 
conditions are symmetric as it explicitly appears in the 
discrete version of the problem (see Ref. [1^). This is 
also emphasized in Fig. yj, where the four matrix elements 
of a classical Green's function Q{t,t') (for i'/T — 0.6) 
are plotted with respect to time t. This example cor- 
responds to a classical orbit of the 2D hydrogen atom 
in a magnetic field having initial and final points on the 
nucleus, namely, the closed orbit having code 0— and 
whose trajectory in (m, v) coordinates is also shown in 
the figure, (see Sec. ^ for all details.) As expected, the 
Green's function vanishes at initial and final times (i.e., 
g{0,t') = Q{T,t') = 0) and for t = t', the derivatives of 
each diagonal elements gii{t',t') (continuous line) and 
G22{t',t') (long dashed line) are discontinuous whereas, 
from property (E5[), the two off-diagonal elements are 
equal (dotted and dashed lines). 



Getting Ci(q, qo,r) by integrating a set of first 
order differential equations 



From Eq. (0), there are three contributions to 
Ci(q, qo,T), namely: 



hiT) = 

I2{T) = 



dtv^'^i{t)g,,{t,t)gki{t,t), 



T pT 



+' T/(3)('+M/(3) (+1 



dtdt' v^'> {t)Vi\ll {t')g,j (t, t)gki {t, t')g„,^ {t', t'), 



10 JO 



+' T/(3)^+\T/(3) W 



(16) 



0"'0 



dtdt' vi;^{t)Vi^Zit')Suit, t')gj„,it, t')gut, t')- 



r 



Even if, in principle, one can compute g{t, t') for any 
(t, t') values using Eq. (p^), direct evaluation of the dou- 
ble integrals /j would be time consuming and numeri- 
cally inefficient using standard integration routines, espe- 
cially because, from its definition, g(t, t') is not a smooth 
function around the line t = t' . In what follows, we will 
show that the preceding integrals can be transformed in 
way such that their values can be obtained integrating a 



set of first order differential equations along the classical 
orbit, in the same way that, for example, the monodromy 
matrix M{T) can be computed. 

Separating t > t' and t < t' contributions in /j , using 
symmetry property (15|) of g{t,t') and that the matrix 
V^^' is fully symmetric under index permutations, one 
gets, after straightforward algebra. 
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Figure 1: Example of a classical Green's function Q{t,t') 
involved in the calculation of the h corrections for the prop- 
agator /^(q, qo,T), for the case q = qo = 0. It is associated 
with the closed orbit 1243 of the 2D hydrogen atom in mag- 
netic field, whose trajectory in (u, v) coordinates is inserted 
in the plot (see Sec. fVllfor all details). This trajectory starts 
and ends at the nucleus, depicted by the black circle. Each 
curve corresponds to a matrix element Qij{t,t') plotted with 
respect to time t, for t' /T — 0.6. As expected from boundary 
conditions (M), the Green's function vanishes at initial and 
final times (i.e. 0(0, t') = g{T,t') = 0) and for t = t' , the 
derivatives of diagonal elements, 6ii(i',i') (continuous line) 
and Q22{t' ,t') (long dashed line), are discontinuous whereas, 
from symmetry property (nsf) (i.e., Q^{t,t') — Q{t',t)), the 
two off-diagonal elements are equal (dotted and dashed lines). 






(17) 



r 



In the preceding expressions the Green's function g{t,t') 
is used only for (i, t') values in the triangle < t' < 
t<T a.nd is formally written g{t, t') = BZ{t)J^{t') (see 
Eq. (p^)), thus separating t and t' contributions: 



This leads us to introduce two intermediate quantities, 
namely, Pp{t) and Qpqr{t) (for p, q and r running from 1 
to/): 



1+ (T) = 2 / dt V^ (t)ay (t, t)B;, it) 
Jo 



+' 1/(3) W 



X / dt'V^ZitVllpitlGmnit'X), 



I^{T)^2 f 
Jo 



dtV^l{t)B.{t)B.{t)B;^[t) 



4-1 T/(3) (,l\ 



(18) 



X / dt' V,'ZitVllp{t')Jlm,{t')Jlnr{t'). 

Jo 



+' 1/(3) (+'\ 



Pp{t)= / dt'ViZ{t')Jnp{t')g„,n{t',t'), 



Qpqr{t) — I dt'Vi^^^{t')Jiip(t')Ji„iq{t')Jinrit') 

Jo 



(19) 



in a way such that I^iT) (and Ii{T)) are solutions of the 
following set of differential equations (besides equations 



for X(i) and M{t)): 



^(4), 



h = vl^i{t)g,,{t,t)gu{t.t), 

Pp = vl±{t)Jnp{t)Q^n{t,t), 



(20) 



Q 



r(3) 



'^lmnWJilpWJir,iq{'t)Jl nr (t) , 



pqr 

I2 = v!^it)B-MB;At)B;,{t)Q.p,r{t) 



with initial conditions /i(0) = ifi^) = -fp(O) = 
QpqriO) = 0. This set of equations, p + 4p + 3f + 3 in 
total (i.e., 33 for a 2D system) is easily integrated using 
any standard method (fourth order Runge-Kutta in the 
present case). As mentioned in the Introduction, the size 
of the preceding differential set is probably not minimal 
and could be reduced by a deeper analysis of the struc- 
ture of these equations. However, it allows a fast and 
easy computation of the correction Ci(q, qo, T): 

• Find a trajectory going from qo to q in time T; 

• Integrate the differential set for X(i) and M{t) 
along the trajectory to obtain the quantity 



• Integrate the set of Eqs. ( pO| ) 
to get the three quantities h, /j 
Ci (q, qo , T) expression. 



along the trajectory 
entering in the 



III. THE GREEN'S FUNCTION G(q,qo,£) 
A. Going from time to energy domain 

since the quantum Green's function G'(q, qo,i?) is re- 
lated to the propagator K{q, qo, T), through a semisided 
Fourier transform, this relation also holds between semi- 
classical contributions arising from each classical orbit, 
more precisely. 



Gi(q,qo,£') 



1 
ih 



+ OC 



dT exp 



^£;TW,(q,qo,r). 



(21) 



Again, a stationary phase approximation is used to per- 
form the integral, which, for a given trajectory going 
from qo to q, selects its total duration To such that 
the classical motion is made at energy E. This opera- 
tion also gives rise to additional terms in h corrections, 
to be summed with Ci (q, qo , T) , and whose explicit ex- 
pressions can be derived starting from Eq. (0) formally 
written as follows [13|: 



Ki{q,(io,T)^ 



1 



■ exp 



^Wi{q,(io,T)-i^J^i 



(27rzft)//2 

+Coiq,qo,T) + thCiiq,qo,T) , (22) 



Co (q, qo , T) being thus the (logarithm of) usual semiclas- 
sical amplitude. Then W^i(q, qo,r) and Co(q, qo,T) are 
systematically expanded around To: 



Wi{q,qo,T) = Wi^'''^ +STWi 



Ai) 



+ 



ST'', 



2 ' 
ST 



"J_W'^ + "4^W^^^\ (23) 
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Co(q,qo,T) = C^"^+<5rC^'^ 



ST^ 



c, 



(2) 



with 6T = [T — To). Terms arising from Ci(T) expan- 
sion would contribute only to h^ correction and can be 
discarded. Performing the imaginary Gaussian integrals 
leads to the additional h corrections: 



Ci^^^(q,qo,ro) 



wl^^c^^ 



1 



2W, 



(2) 



K0'+^' 



w, 



(4) 



wl 



(3) 



{Wl-^)' 8(Tyf))^ ''(Wl-^f 



(24) 



The preceding formula is similar to the one in 
Ref. [n3 , where the authors have expressed the coefficient 
Ci^^q, qo. To) in terms of derivatives of amplitude and 
action with respect to energy E. The full expression of 
G;(q, qo, E) is then given by 



G((q,qo,-B) 



27r 



1 



2T:ih) 
X exp 



(/+l)/2 



VF/'MetJi(To) 



1/2 



-Si{(\,qa,E)~i-vi 



{l + ift[Ci(q,qo,To)-HGf^^(q,qo,To)] -1-0(^2)1^ 



where Si (q, qo , E) is the reduced action and 

.(2) 



i^i = n 



if wl ' > 0, 

(2) 



Di^ui + 1 if wr> < 0. 



(25) 



(26) 



B. Getting C]^~*^(q, qo,To) by integrating a set of 
first order differential equations 



In Sec. [IC, we have shown that Gi(q, qo,To) can be 
computed by integrating a set of differential equations 
along the classical orbit going from qo to q in time 
To. In this section we will show that it is also true 
for G^^^(q, qo. To), which involves derivatives of both 
Wi{q,qo,T) and det Ji(T) with respect to T. 

For all T, we have the following functional relation (qo 
and q being fixed): 



dWiiq,qo,T) 
dT 



= -S(q,qo,T), 



(27) 



where E{q, qo, T) is the energy of the classical trajectory, 
q(i, T), going from qo to q in time T, that is, the value of 
the Hamiltonian H taken at any point on the correspond- 
ing phase space trajectory X(i, T) — (q(i, T), p{t, T)) . 

Writing T — Tq + 6T, the Taylor expansion of 
iJ(X(i,To + (5T)) is easily deduced from the Taylor ex- 
pansion of X(t, To + ST) around the reference trajectory 
X(t,To) (noted hereafter as X(")(t))): 



where again all derivatives of H are evaluated at X(°^ (t) . 
Thus, we are facing three differential sets of the form 
X'*) — Sff^^^X'*^ -I- SY^*) (i.e,. nonhomogeneous hnear 
differential equations), with the important property that 
the vector Y^'^ only depends on vectors X^-') with j < i, 
so that they can be solved one after the other. Solutions 
of these nonhomogeneous linear differential equations are 
expressed with the monodromy matrix M^*^': 



X(i, To + 5T) = X^°'> (t) + ST X^i) (i) 



+ ^x(^)(t) 



6 



X(3)(i) 



(28) 



and from which one obtains the higher derivatives of the 
classical action Wi at T = Tq: 



X(i)(t) =A'f(")(i)X(i)(0,) 

X(2)(t)=Af(°)(i)X(2)(0)+F(2)(t), 

X(3)(i)=Af(")(i)X(3)(0)-KF(3)(t). 



(31) 



W^/^) = -if(x(o)(i)), 



^(4) ^ _ /^(3)^(1) 



-Aj A^. Aj, U-.^. 



Among the 3 x (2/) dimensional space of solutions given 
(29) by preceding expressions, the relevant one is selected by 
transposing on initial values X*^*)(0) (for i = 1,2,3) the 
two boundary conditions 



where all derivatives of H are evaluated at X*^"^(t). 

Equations for X("^(i) are deduced from Hamilton's 
equations governing X(t, T) evolution: 



q(0,To + 5T) = qo and q{To + ST,n + ST) 



q 

(32) 



X 



(1) 



(2) v(l) 






^i - ^ij^jk ^k 

X^'^ = ^^jH^Xl^ 
t '■J jk k 



V H-(3) T^(l) y(l) 



(30) 



r(4) 



(i)y(i)y(i) 






Introducing position q'*^ and momentum p*^*) parts for 
vectors X'^' , the Taylor expansion of the preceding equa- 
tions leads to the following boundary conditions: 



f q(i)(0) = 

q(2)(0) = 

lq('nO) = 



and 



qW(ro) = -qW(To), 

q(2)(ro) = -qW(To)-2q(i)(To), 

q(^Hro) = -q^"HTo) - 3q(i)(To) - 3q(2)(To). 



(33) 



r 



Thus, the initial values p*-*-'(0) are implicitly determined 
by the final values q^*^(To), through the integral expres- 
y), which for X^^^ reads 



sions 






J2{To) Ji(To) 
Mn) Ji(To) 



p(i)(0) 



(34) 



showing thus that p^^^O) = -Ji\To)q^°'>{To). 

Then F'^^^Tq) and F^^'^{Tq) are easily computed by 
integrating sets of differential equations obtained from 
Eq. (SO), allowing us to derive p''^-'(0) and p('^)(0) values 



from Eq. (pi|), solving systems similar to Eq. (p4|): 

' p(2)(0) = -J^\To) (q(")(To) + 2q(i)(ro) +f(2)(To) 
p(3)(0) = -Ji\To) (q^'^HTo) + 3q(i)(To) 

+3q(2)(To)+f(3)(ro)V 



(35) 

where we have introduced the notation (f^*-*, g^*-*) for vec- 
tors fW. Quantities like qf^H^o), q^^H^o), and q^^H^o) 



can also be expressed in terms of X^"^ (Tq) and its deriva- 
tives. 

At this point, from the values of the three vectors 
X(*^(To) and using Eqs. ( |29| ) at time Tq, all derivatives 
T4^("^ of the classical action can be computed. 

We now explain how to compute derivatives of 
det Ji (T) . More precisely one has to calculate the 

two coefficients Cq and Cq , which are derivatives of 
— In -^/l det Ji(r)|, so that, using the well-known formula 



A(l„|de.7|)^Tr(j-.g 



(36) 



(J being any (invertible) matrix), expressions of Cq and 



(1) 



Cq become 



C, 



(1) 



-itr( Jr^(To 



dT 



Co --2^1 Ki ™^T^ 



-J.-\n)'-^j^\n)''^^^^^ 



dT 



dT 



(37) 



where dJi(To)/dT means derivative of Ji(ro) when 
changing total time T (and thus the classical orbit), 
which must not be confused with Ji (time derivative 
of Ji along a given classical orbit). Ji{T) being the 
(/ X /) upper right submatrix of the monodromy matrix 
AI{T), d" Ji(ro)/(iT" is also stored at the same position 
in matrix d" M (Tq) / dT^^ , for which we will derive gen- 
eral expressions. For this purpose, we first introduce the 
explicit notation M(t,T), representing the value of the 
monodromy matrix at time t along the orbit going from 
qo to q in time T. Writing T = Tq + 6T, the Taylor 
expansion oi M[t, T) for a given time t reads 



M{t,To + 5T) = M^^^\t) + 5T M^^\t) + — -M(2)(t), 



(38) 



where IvT^'it) is the monodromy matrix along the ref- 
erence orbit (i.e., going from qg to q in time Tq). Then 
dM{To)/dT and d'^M{To)/dT^ arc the Taylor coefBcients 
of monodromy matrix M{Tq + 5T, Tq + ST) and thus have 
the following expression: 



^^=Af(°)(ro) + MW(To), 



dT 



f^^ ^ Af (o)(ro) + 2M(i)(ro) + Af(^)(To). 



(39) 



Equations governing Af ('^(t) evolution are easily deduced 



from the one for M{t, T): 



Mi;' = E,fc 



H^M^+H^Xl^M^ 



M. 



(2) 



i^lf<' 



^HJ^t^^Mlp 






(40) 

with initial conditions Af(i)(0) = A/(2)(0) = 0. Obvi- 
ously these equations are similar to those governing X^*) 
evolution, so that M'^^^Tq) and M'-^^Tq) values will be 
obtained by integrating similar differential sets. Actu- 
ally, it can be shown that all these sets (for both X*^*' 
and Af'^*)) can be concatened in only one (larger) set of 
differential equations, whose integration can be done at 
once. 

Finally, gathering all quantities in Eq. (p9), the two 
matrices dJi{To)/dT and d^ Ji{Tq) / dT'^ are inserted in 
Eq. (g^ thus giving values for Cq and Cq , which, 

in) 

along with the values for Wj , allow us to compute the 
numerical value for Cf ^^(q, qo,TQ). 

Obviously, the number of equations in the preceding 
differential sets can be reduced, especially for Hamil- 
tonian separating into kinetic and potential energy, for 
which -ff,L; and ^L;„ coefficients are nonvanishing only 
when 1 < j,k,l,m < f. However, these sets are straight- 
forward to implement and need only a small amount of 
CPU time to be solved using any conventional integrator 
(4th order Runge-Kutta in the present case). 



IV. TRACE OF THE PROPAGATOR K{T) 

The diagonal elements -ftr(qo, qo, T) of the propagator 
are related to classical orbits starting from qo and return- 
ing to this point after time T, i.e., closed orbits. Sum- 
ming all these diagonal elements, that is performing the 
integral / dqo K{qo,qfj,T), will select, through another 
stationary phase approximation, closed orbits for which 
initial and final momentum are equal: periodic orbits, h 
corrections to leading order of the semiclassical contri- 
bution to K{T) from each periodic orbit can be derived 
following the same scheme, previously used for the prop- 
agator itself ||, |l|| . 



A. Feynman path integral 

Adding the integral over the initial and final positions 
in Eq. (|) yields H 



K{T) = / dqo dqi dq2 , • ' ' , dqN-ii2TTihAt) 



X exp 






At 



-Nf/2 



q„ At-t-C'(At) 



(41) 
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with q^ = qo. 

The stationary phase approximation around a given 
periodic orbit (lf{t) is made exphcit when replacing the 
preceding Nf integral with ||l^ 



ydgljd|o^d|id6,---,d|iV-i 



(42) 



with ^n = Qn — qf{nAt). For n = (i.e., initial posi- 
tion), only deviations perpendicular to the periodic or- 
bit ^Q have been introduced because the classical action 
Wi(qo, qo, T) is constant along the orbit (depicted by gjj). 
The contribution Ki{T) of this periodic orbit to K{T) 
then reads [Ol 



Ki{T) = 



N 



Nf/2 



2TrihT^ 
X d|2 , • • 



exp ( -Wi 



dglj d^o"- d|i 



,d|Ar_i exp ( —W^ab£,a^b 



-W^abc^a^b^c 



72ft2 



W,abcW,def£,a^b^c^d£,e£,f 



W^abcdS,a^b^c^d 



(43) 



where ^a — ^oi when a — (0, i) and ^a = when a = 
(0,0). Wi is, in the large N limit, the classical action of 
the periodic orbit. Full expressions for W^ab, W.abc, and 
W^abcd can be found in Ref. Q . 

Then, the next step would consist of performing all 
imaginary Gaussian integrals, leaving out the integral 
along the orbit. However, in the preceding coordinate 
transformation (|4^), there is an hidden subtlety, affect- 
ing only h corrections, which explains why probably it is 
not mentioned in usual textbooks ||l|, g], where authors 
are only looking at leading semiclassical amplitudes. 

Actually, the problem is that the integral over q^ cor- 
responds to the length of the classical orbit, only when 
^(j- = 0; for a nonzero value, it will correspond to in- 
tegration on a closed curve, slightly displaced from the 



original trajectory, whose length will thus depend on the 
$,Q value. To enlighten this, let us suppose that we 
have a bidimensional system, for which one periodic or- 
bit is a circle of radius i?o, traveled at constant speed 
Vo ~ 2'kRq/T. The coordinate transformation is then 
easily made using polar coordinates (r, 9): 



r = Ro- ^0 



(44) 



The negative sign appears to preserve orientation. The 
volume element dx dy becomes 



dx dy = rde dr = (i?o - ^o)d0 d^^, 



(45) 



which shows that, in this case, dqQ is not simply RodO, 
the length on the periodic orbit, but is given by 

dql = (i?o - io)dO + Rode. (46) 

This simple example shows actually that the variable ^q 
is not independent of ^q, whereas is. 

For a general system, the variable that can play the 9 
role is actually the time f , whose variation domain [0, T] 
is fixed and then obviously independent of $,q . Thus one 

has to generalize the relation dq'^ = jq'^'ldio, valid only 
on the periodic orbit. This is done by writing explicitly 
the coordinate transformation q ^ (^o, ^o" )■ 



q^'(to) + eo>.(io), 



(47) 



where ni{to) are / — 1 orthogonal unit vectors lying in 
the plane perpendicular to the periodic orbit at time tg. 
The Jacobian of the transformation reads 



det 



9q 



9(to,IJ 



det [cf^ + ^(ji n,, ni, • • • , n/_i] 



1^0 ■ q' 



(48) 



Inserting the volume element in Eq. (^3|) , the contribu- 
tion Ki(T) of the periodic orbit now reads, keeping only 
terms giving rise to h corrections. 



KiiT) = 



N 



27rihT 



Nf/2 



e^p['-Wij \cf'\dto d|o^ d|i d^2 , • • • ,d^N-i exp (^W^abUb 



\ , 0"Krf * w C C C . * W C C C C . ^ VjW,abc^^^aib£.c 

1 + 1^^ + -^W .abcS.aS.b£.c + -;^W ^abcdUKbUU + ^- 



1 



|qcl|2 



72ft2 



W^abcW.defiaibicidieif 



, (49) 



r 



where we have seen that q'^' = ^dqV and we have intro- are extracted when comparing the detailed expression of 



duced the index d for (0, j). 

As explained in Ref. |l3], the imaginary Gaussian in- 
tegrals can be expressed in terms of another classical 
Green's functions Q(t,t'), whose boundary conditions 



W,ab with the discrete version of the Jacobi-Hill operator 
V, see Eq. (^. Especially, it can be shown that, in the 



large N limit, they become 

Vt,g{o,t')^Vt„giT,t') = o, 



Vt'e[0,T], (50) 



where we have introduced Viq the projector along the 
periodic orbit at time to and Qt^ = 1 — Vtg ■ In Ref. Q , 
only the f^ + f boundary conditions corresponding to 
the first two lines were given, whereas the P — f ones 
corresponding to the last line were missing. 

Gathering all results and taking the large N limit in 
Eq. (MS) , the contribution of the given periodic orbit to 
the trace of the propagator reads as follows: 



Ki{T)^Kf'\T){l + ih 



dtoCi{T,to)+0{h^ 



(51) 



K, JT) being the usual semiclassical leading order R, 



k['\t) 



1 



T 



\/2^ I^ETdet (m(r) - 1)1^/2 
X exp 



'jWi{T) - i^tii + isgn BeT 



(52) 



where Wi (T) is the classical action of the periodic orbit 
and /i; its Maslov index. 

The first h correction Ci(T) to kI°\t) is then ob- 
tained by averaging over the time to (i.e. over the full 
periodic orbit) the coefhcient Ci{T,to), given by 



Ci (T, to) = lJ^ dt v^i {t)g,, (i, t)gki {t, t) + i ^^^j^ ^ dt v^ {t)gik (o, t)g,, (t, t) 



T pT 



+ ^J J dtdt' l/>-t^(t)T^J:i(iO[3ey(i,t)eM(t,i0^mn(i',t') + 2a.i(i,t')ejm(i,i')efcn(i,i')], (53) 



r 



where to represents thus the position qo on the periodic 
orbit at which boundary conditions (pW on the classical 
Green's function g{t,t') are applied, qo is also the ini- 
tial (and final) position on the periodic orbit for classical 
motions corresponding to times t and t' entered in the 
preceding expression. 



B. Classical Green's function 



As in Sec. |IIB| , where expressions for classical Green's 
functions for the propagator K{c[,qo,T) where derived, 
we introduce the g±{t,t') notations and A±{t'), B±{t') 
matrices. Using all boundary conditions (at times t = t' , 
t — and t — T) gives rise to the following equation: 



"1 

Qto 


(Af (T) ~ t2f) ( 






■ 1 " 






[O Qt„\ 



A At') 
B-{t') 

M{T) 



-Jlit') 
Jlit') 



(54) 



The preceding set of linear equations, formally written 
AX = B, cannot be solved directly because the (2/ x 2/) 
matrix A is obviously singular. More precisely, existence 
and number of solutions for the system yl x = b are de- 
termined by the two following properties: 



• Solutions exists if for all vectors y such that A y = 
0, then y • b = . 

• If the preceding condition is fulfilled, and if x is a 
solution, then for all vectors xq such that Axo =0, 
x-|-Xo is also a solution, showing that the dimension 
of the solution space is that of the nullspace of A. 

In the present case, equation A^ y = leads to, cither 



1 

Qto 



y oc 





q(io) 



(55) 



or 



(Af (T) - l2/)y = with y = S 



1 

2*0 



y^O. (56) 



For a generic unstable periodic orbit, the eigenspace as- 
sociated with the eigenvalue 1 of M{T) (T being the pe- 
riod), is of dimension one and is spanned by the vec- 
tor parallel to the flow (see the Appendix) X(io) = 
(q(to),p(^o)), so that, in the second case, one gets y ex 
X(io) and y is a solution of: 



1 

Sto 



y ex 



-P(io) 
q(io) 



(57) 



which is impossible unless q(to) = 0, which, for Hamil- 
tonian separating into kinetic and potential energies. 
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corresponds to a self-retracing periodic orbit, for which 
a sUghtly modified approach should be developed [[l8[ . 
Nevertheless, this case is peculiar, and we will suppose 
in the rest of the section that q never vanishes along the 
periodic orbit in consideration. 

Thus, the nuUspace of A^ being one-dimensional and 
spanned by the vector (0,q(to)), Eq. (M) immediately 
shows that for any column of matrix B, we get (0, q(io)) • 
Bi — 0, thus fulfilling the first condition. Denoting Xq 
as a solution of Eq. (p4|), which can be easily obtained 
using singular value decomposition (SVD) of matrix A, 
and the nullspace of M(T) — 1 being spanned by X(fo), 
the general solution of Eq. (M) reads 



A. — /Xn 



aiX(io), a2X(to), • • • , a/X(to) 



(58) 



where a^ are unknown real parameters still to be deter- 
mined. Actually, in Eq. (|5J) one boundary condition has 
not been taken into account, namely that VtoG-{0, t') — 
which, using that the projector Vta reads 



(^to) 



q(to)q^(to) 

|q(io)P 






(59) 



allows us to get ai values and, from that, the final ex- 
pression 



A-{t') 
B-{t') 



— Xq 



1 



|q(io) 



q(io)q^ 
p(to)q^ 



(io) 
{to) 



Xn 



(60) 



which, of course, is now independent of the particular 
solution Xq. 

Whereas in the case of the propagator K{q, qg, T), for 



which we were able to give an explicit expression (14), 
the classical Green's function associated with the trace 
of the propagator K{T) is only defined trough a linear 
system (pj), which nevertheless allows us to obtain its 
numerical value for any (i, i'). Although it clearly ap- 
pears that matrix D expression (see Ref. [O ) is symmet- 
ric, meaning that the classical Green's function must ful- 
fill the property Q^{t,t') = Q{t',t), getting the later di- 
rectly from Eq. ( p4|) is not obvious. However, in the case 
of the 2D hydrogen in a magnetic field (see Sec. W^ for 
all details), we have numerically checked that the prop- 
erty holds. For example, in Fig. g the four coefficients of 
classical Green's function G{t,t') (for t' /T = 0.3) of the 
periodic orbit 1234 are plotted with respect to time t. 
The starting point to on the periodic orbit is depicted by 
the cross. Actually, we have plotted the coefficient of the 
rotated matrix Q{t, i'), such that its first row corresponds 
to the direction parallel to the orbit; Gii{t,t') (continu- 
ous line) and Gi2it, t') (dotted fine) are thus equal to zero 
for initial [t = 0) and final (i = T) points. The other 
boundary conditions can also be verified in the figure: 
the dashed fine {G2i{t,t')) (respectively the long dashed 
line, G22{t,t')) has not only the same value at initial and 
final time, but also the same slope, which means that 



o" 




-0.5 



Figure 2: Example of a classical Green's function Q(t,t') in- 
volved in the calculation of the h corrections for the trace of 
the propagator K{T). It is associated with the periodic or- 
bit 1234 of the 2D hydrogen atom in a magnetic field, whose 
trajectory in (it, v) coordinates is inserted in the plot (see 
Sec. Knjfor all details). The black circle depicts the nucleus, 
whereas the cross corresponds to the initial and final points 
on the periodic orbit at which Q{t,t') fulfills the boundary 
conditions (pOl). Each curve corresponds to a matrix element 
Gij{t,t') plotted with respect to time t, for t' /T = 0.3. Ac- 
tually, we have plotted the coefficient of the rotated matrix 
Q{t,t'), such that its first row corresponds to the direction 
parallel to the orbit; Q\i{t,t') (continuous line) and Qi2{t,t') 
(dotted line) are thus equal to zero for initial (f = 0) and final 
(i = T) points. The other boundary conditions can also be 
verified in the figure: the dashed line (C/2i(i, i')) (respectively 
the long dashed line, G22it, t')) has not only the same value at 
initial and final time, but also the same slope, which means 
that Q2i(t,t') (respectively. G22{t,t')) and its time deriva- 
tive fulfills the periodic boundary conditions (BQ). Finally, 
for t = t' , the off-diagonal coefficients Gi2{t',t') (dotted line) 
and G2i{t' ,t') (dashed line) are equal, as expected from the 
symmetry property G^{t,t') = Q{t' ,t) . 



G2i{t,t') (respectively G22{t,t')) and its time derivative 
fulfill the periodic boundary conditions (50). Finally, for 
t — t' , the off-diagonal coefficients Gi2{t' , t') (dotted line) 
and G2i{t',t') (dashed fine) are equal, as expected from 
the symmetry property. 



Getting Ci{T,to) by integrating a set of first 
order differential equations 



As seen previously (see Sec. 11 C), we will explain how 
the numerical value of coefficients Ci{T,to) can be ob- 
tained by integrating a set of differential equation, using 
the standard Runge-Kutta method. There are now four 
contributions to Ci(T, io), namely 
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.(4) 



/i(T)= / dtv^;^i{t)g,,{t,t)gkiit,t) 



r{3) 



ii{T)= / dtv^'^^{t)gik{Q,t)g,j{t,t) 



liiT)^ fj^tdt' Vifl{t)vtt(t') 
JaJQ 

xga{t,t')gj„,{t,t')gkn{t,t'). 



(61) 



The two main difficulties now are that g{t,t') does not factorize anymore in a product of matrix at time t and a 
matrix at time t' , nor does the symmetric property g^ {t, t') = g{t\ t) explicitly appear (even if we have numerically 
checked that it is fulfilled). Nevertheless, as seen previously, separating {t > t') and (i < t') contributions in /2 (^) 

expressions and introducing four quantities Pp for 1 < i < 4, allows us to compute I2 ^ by integrating the following 
set of differential equations from t = to T (besides equations for X(t) and M{t)): 



p(i) 
p 

p(2) 
V 

it 






p(3) 
P 

p(4) 
P 



Jup{t)Vi^l{t)g^^{t,t), 

J2lp{t)V}:±{t)gran{t,t), 



vS it)g., {t, t)A;, {t)p^^) (t) + v^ (t)g,, (t, t)B;, (t)pw (i) 



(62) 



with vanishing initial conditions for Pp and /j . For each time step, one must compute matrices A- and B- (and 
from there matrices ^4+ and i?+), solving the linear system described in the previous section, using singular value 
decomposition of matrix A, which, being independent of t, is done before starting the Runge-Kutta integration. 
Skipping intermediate steps, the differential equations leading to /;f (T) computation reads as follow, introducing 

(i) 

another eight quantities Qpqr'. 

' Q'p'Jr^V,'lit)A+{t)A+^{t)AXM. 



^pqr 
^pqr 

of-*) 



Qpqr = ^lmnWJ^lpi^)'^^rnq{t)J2„r{'t), 
Qpqr = VlmnWJ'^lpWJ2mq{t)Jlnr{t), 
Qpqr = Vlmni^)J2lpit)Jimq{t)Jlnr{t), 



Qpqr — ^lmnV')dllpV-)'JlmqV-)'Jlnr 



(t), 



Vi^^l{t)A+{t)B+^{t)B+M, 
Vi^llit)B+{t)B+^{t)B+^{t), 

= VS{t)J2.p{t)J2,q{t)J2kr{t)Q^pqi{t) + 3V!^{t)J2,p{t)J2,q{t)J,Ut)Ql?ql{t) 
+ SVg!{t)J2rp{t)J^,q{t)J,Ut)QilUt) + Vg!{t)Jup{t)Jl,q{t)J^Ut)QiU*) 

+ vSit)A;MA-,it)A;,it)Q(^\lit) + SKf^it)A;MA-^it)B;,it)QllUt) 
+ ^V^^it)Apdt)B;,{t)B;,{t)Q(;j^{t) + vl^{t)B;^{t)B-^{t)B;,{t)Qll{t) 

I 



(63) 



with vanishing initial conditions for Qpqr and I2 Finally, namely 
one must add equations leading to /; and /i computation, 



A = ^S(t)&,(i,t)afci(t,i), 
ii = vSit)Ai,{t)g.,{t,t), 



(64) 
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where we have used Gik{0,t) = ^^.{t). Taking into ac- 
count equations for X(i) and M{t), this gives rise to a 
total of 8p + 4/2 + 7/ + 3 equations, that is 97 for a 2D 
system. 

In practice, having found a periodic orbit and for a 
given to along this orbit, the coefficient Ci{T,to) is com- 
puted in two steps: 

• Achieve the SVD decomposition of the matrix A, 
appearing on the left-hand side of Eq. (|5J), and 
compute the projector matrix appearing on the 
right-hand side of Eq. (JGy) . 



• Integrate the differential set (|63|) along the peri- 
odic orbit (starting at point depicted by to). At 
any time t, use the preceding SVD decomposition 
to obtain a solution Aq and the projector matrix 
to get the true solution {A_{t),B_{t)) and thus 
(A+(i),S+(i)), using Eq. @. 

Finally, the coefficient Ci{T,to), being a smooth func- 
tion of to, the average over time to, leading to the h cor- 
rection term Ci (T) , can be handled by any conventional 
integrator. 



V. TRACE OF THE GREEN'S FUNCTION G{E) 

Steps leading to the semiclassical contribution Gi{E) 
from a given periodic orbit to the trace of the 
Green's function G{E) are identical to those giving the 
G;(q, qo,£') expression, so that Gi{E) reads 



Gi{E) 



■ exp 



-ME) 






«^|det(m(ro)-l)|^/^ 
X {l + ift [Ci(To) -f Cf ^^(To)] + Oih")], (65) 

where G'[~*^{T) is given by 



Cr^(To) 



2W, 



(3)r^(i) 



(2) 



KO +^0 



(2) 



wr'G, 



w, 



(4) 



(i^f)) s{w! 



(2) 






(66) 



Wi (respectively Cg ) are the Taylor coefficients of the 
Wi{T) (resppectively Go{T)) expansion around Tq. 

Computation of W[ is much the same as in the 
Green's function case, because the functional relation 

dWi{T) 



dT 



-E{T) 



(67) 



still holds for a given periodic orbit, E{T) being its en- 
ergy as function of its period, which is still given by 
the value of the Hamiltonian H taken at any point 
on the corresponding phase space trajectory X(t, T) = 



'^"^ coefficients (Eq. ^ and for X(") (t) 
The only differences with the preced- 



(q(t, T), p(i, T)). Thus, the Taylor expansion of X(t, T) 
around the periodic orbit X(i, To), will lead to the same 

expressions for Wi 
equations (Eq. |3^) . 
ing section arise from the boundary conditions fulfilled by 
X(")(i), deduced from the equation X(0,T) = X{T,T), 
i.e., X(i,T) is a periodic orbit of period T. The Taylor 
expansion of this relation leads to the following condi- 
tions: 

x(i)(o) = x(i)(To) + x(")(ro), 

X(2) (0) = X(2) (To) + X(") (To) + 2X(i) (To) , 

X(3) (0) = X(3) (To) + x'"^ (To) + 3X(i) (To) + SX^^) (To) 

(68) 



Solutions of the differential set (pG) still have the fol- 
lowing formal expressions (p^), which, inserted in the 
boundary conditions (|68|), leads to equations on X(*)(0) 
only 

'(l-Af(To))X(i)(0) = xW(To), 

(1 - Af(To))x(2)(0) = X(")(To) + 2±'-'\To) + f'-^\To), 

(1 - Af(To))X(3)(0) = X^"'(To) +3X(i)(To) 

+ 3X(2)(To)+F(3)(ro). 

(69) 

The matrix 1 — M(To) being singular, solving the preced- 
ing linear equations need additional discussion, which, for 
simplicity, will focus on X(^)(0) only. First, the nullspace 
of 1 - M(To)T is spanned by SX(°)(To), which is obvi- 
ously orthogonal to X'^'^^(ro), the right hand-side of the 
equation for X^^)(0), thus showing that this equation ad- 
mits solutions. Then, the nullspace of 1 — M(ro) being 
spanned by X(°)(ro), the whole set of solutions reads 



X(i)(0)-xW(0) + aXW(ro), 



(70) 



where Xq (0) is a particular solution of the equation. 
Actually, the term aX(°)(ro) corresponds to an displace- 
ment of the initial conditions along the flow, which, of 
course, gives back the same periodic orbit (at first order 
in T— To). We thus expect that this term has a vanishing 

(2) 

contribution to W^ , which is easily verified when insert- 

ing the general solution in the Wi ' expression (taken at 
time t = To): 



wl^ = 



(x«(0)-XW(To) 

+aX(")(To)) •VHfx(")(To: 
= -X[,')(0)-Vi/(x(")(To)' 



(71) 



because of the Hamilton's equations X^*'^(ro) 

i]Vff(x(o)(ro)). 
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These two properties also hold in the cases of X^^^ (0) 
and X(^^(0), but are slightly more complicated to estab- 
lish because the right-hand sides of the equations involve 
F(*)(To) and derivatives of X(^)(ro). 

Thus, integrating the same differential sets that were 
used for G(q, qo, E), one is able to compute the first four 

derivatives of the action, Wi , with respect to the period. 

Starting from the Cg (T) expression 



which means that one is left with the calculation of 
derivatives of In | det (m(T) — 1)| with respect to the pe- 
riod T. As shown in the Appendix, det {m{T) — 1) is 
given by the determinant of the 2/ x 2/ matrix N{T) 
defined as follows: 



N{T) = M{T) - (1 - V\\ (T) - V^{T)) , 



(74) 



cI^\t) = In T - - In \dET\ - - In | det {m{T) 1)| 

(72) 

and using the fact that OeT — l/drE = —1/d^Wi, one 
obtains 



L-O (Joj-TTT + o 



Cf(To) 



(3) 



To 2 v^^(2) 2 dT 



l^ln|det(m(r)-l)|. 



IW, 



(4) 



1 w, 



(3)' 






(2) 



2w^ ' '^ \Wi 
i^ln|det(m(T)-l)|, 



(73) 



where we have introduced P\\{T) (respectively V±{T)) 
the projector on the direction parallel to the flow (respec- 
tively perpendicular to the energy shell), more precisely, 
the 7^11 (T) and V± (T) expressions are 



Vn = en 



and V± = e± ■ el = -ST'iii:, 



(75) 



where eii is the unit vector tangent to the flow at initial 
(and thus final) time and e^ = Se||. Now, using again 
formula (p6h, derivatives of det {m{T) — 1) with respect 
to the period read 



A(det(m(T)-ll)) = tr(A^(ro)-^^ 

d^ 
dT^ 



(det {m{T) - 1)) = tr (n-\To) 



dT^ 



iV(T„)-^^A.(To)-^^(^°^ 



dT 



dT 



(76) 



with 



dN{To) _ dMiTo) , dVniTo) 



dT 



dT 



dT^ 



dT^ 



dT 



d^N{To) _dHl{To) (PV\\{To) 



dT^ 



dT 
^ dT^ 



(77) 



S. 



As seen previously (Sec. IIIB), 



dMjTo) 
dT 



and — ^j.'i ° ' are expressed in terms of the coefficients M^'^'{t) of the Taylor 
expansion of the monodromy matrix M{t,T) (associated with the periodic orbit X(t,T) of period T) around the 
periodic orbit X(°)(i) of period To, see Eq. (ph. 



Inserting the Taylor expansion of X(T) around To in the V\\ (T) expression, namely. 



'P\\iT) 



|X(T)|| 



-X(T) • X(T)^ 



(78) 
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one obtains the derivatives of V\\ (T) with respect to T: 

( dvwin) _ 1 



dT 



|X(")||2 



±w .±(oV +±io) .±(1)^) _2^^-^^^,rn{To) 

J ||X(")||2 "^ ' 



fZ01 ^ _^±_ (x(2) . x(o)" + X(o) . X(2)" + 2X(i) . X(i)' 

dr2 ||x(o)||2 V 



/ (±io)\±wY x(")^.x(2) x(i)^.x(i)^ 

8-!^ ■■ — ^ 2 .- 2 ^ — ^ 



|X(")||4 



V 

l|X(0)||4 



|X(o)||2 



||X(")||2 



(x(i).XW"+X(o).X(i)"), 



(79) 



Vn(To) 



where all X^*^ are evaluated at time t = 0. 

Gathering the preceding expressions into Eq. d7q) al- 
lows us to compute Indet (m{T) — 1) derivatives, which, 
inserted together with derivatives of the action, in 
Eq. d73) gives the numerical values for Cq (Tq) and 

(2) 

Cq {Tq), which finally leads to the additional h correc- 
tion Cf^^iTo). 



A. Quantum and classical properties 

In atomic units the Hamiltonian of the 2D hydrogen 
in magnetic field reads 



.^ip= 



1 



1 



\/x^ + y'^ 



2 2 

7 y , 



(80) 



where 7 = B/Bq, with Bq = 2.35 x lO^T. The classical 
counterpart of this Hamiltonian has a scaling property, 
that is, if we define new variables by 



(81) 




VI. APPLICATION TO THE 2D HYDROGEN 
ATOM IN A MAGNETIC FIELD 



The hydrogen atom is one example of quantum sys- 
tem whose classical counterpart depicts a chaotic behav- 
ior and has been widely studied (see for, e.g., Ref. Q] 
for a complete review). It has now become a very useful 
tool for testing new ideas and tools in the quantum chaos 
area, both on the semiclassical pG, E3] or universality |24| 
points of view, especially because computing very highly 
excited states has become a standard task on a regu- 
lar workstation, allowing the semiclassical regime to be 
reached easily. Even if one would have preferred to work 
with the real hydrogen atom (i.e., the three-dimensional 
one), in this paper we will focus on the two dimensional 
hydrogen atom in a magnetic field, because taking into 
account invariance by rotation around the magnetic field, 
gives rise to centrifugal terms in the Hamiltonian (typi- 
cally L2fi,2/2r2) which would also contribute to h correc- 
tions and would need a study on its own. One must also 
notice that, even if the classical dynamics are identical 
for both cases, the fact that the magnetic field axis is no 
longer a rotation axis in the 2D case gives rise to slight 
modifications in the Maslov indices [Esl ^ ^. 



we obtain a new Hamiltonian H given by 



H = j-^'-'H = 



2 y^x^ + y 



r,2 



(82) 



which docs not depend on 7 anymore. The classical dy- 
namics of this Hamiltonian is entirely fixed by the scaled 
energy e given by: 



- j-^/^E. 



e==7 



(83) 



All properties of the classical trajectories of the original 
Hamiltonian can be deduced from the scaled dynamics 
using the scaling transformation ( |8l| ) . From the quantum 
point of view, this scaling introduces an effective h value, 
which is easily seen on the scaled Schrodinger equation, 
Htjj — e%p, for a fixed scaled energy e: 



1' 



2/3 



-Af- 



1 



y' 



y/x^ + ]p 



ijj — eip. 



(84) 



Thus, the effective h is given by 7^/'^ and so at a fixed 
value of the scaled energy e, the semiclassical limit is 
obtained when 7 tends to 0. 

The singularity in the classical equations of motion due 
to the divergence of the Coulomb potential at r = 
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is regularized using the semiparabolic coordinates (u — 
\/f + a;, V = \Jr — x) , giving rise to the following effective 
classical Hamiltonian 
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provided it is far enough from any bifurcation and that 
q and qo are not conjugate points for this trajectory. 
Amplitudes and phases being defined by 



n=\pl + \pI - e{u' + v') + lu\\u' + v% (85) 

the trajectories corresponding to the original problem are 
obtained when fixing total energy TC = 2. The associated 
quantum Hamiltonian reads 



A, 



1 



A2 



Wl'> det Ji(ro) 



1/2 ■ 



H{h) = 



d^ 



92 



2 V du'^ dv-^ 



1 



/ Z , Z\ t Z Z/ Z , Z\ 

(86) 



Si = S'(q,qo,A), 



(91) 



which separates into kinetic and potential energy, so that 
the semiclassical formula derived in the preceding sec- 
tions applied to the associated quantum Green's function 
G{z, h), the hydrogen in a magnetic field being recovered 
for z = 2 (actually z/2 corresponds to the nucleus charge) 



G{z,h) 



n{h) 



E 



\r,h){TM 

Z — Xrih) ' 



(87) 



I C, = Ci(q,qo,ro) + Cf-'^(q,qo,To). 

Neglecting h corrections in Eq. (|9y) , the Fourier trans- 
form with respect to the variable C = 1/^ of the following 
function: 



('27r^3/2 1 

9oi0 - ^-^ X --Im(q|G(A,C)|qo) = 

^^^ Y, AA^)^r.cMC'^'S (A - A,(C)) (92) 



where |r, H) is an (normalized) eigenvector of 'H{fi) for 
the eigenenergy Ai-(fi), t representing the set of quan- 
tum labels, i.e., level number and symmetry proper- 
ties (see below), describing |t, ?i). The matrix element 
{q\G{z, ?i)|qo), where q = {u, v) then reads 



(q|G(z,n)|qo) =Y^rA'i)'^r,h{^a) 



1 



z — \t{K) ' 



where V'T,?i(q) = (qI''', ft) has been supposed to be real, 
with T-LiK) being invariant under p — > —p. Taking z — \ 
on the real axis, the imaginary part of (q|G(z, ?i)|qo), 
becomes 

--Im (q|G(A, ;i)|qo) = V V'r,n(q)V'r,;.(qo)'5 (A - K{h)) 

(89) 

to which any classical path going from q to qp at energy 
A, gives the following contribution (see Eq. (pa)): 



ilm(q|G(A,;i)|qo)j 

IT 



^^Jp7j^,|cosQ5, + 0, 



-hCi sin ( -Si + (pi 



(90) 



will depict peaks at the classical actions Si/2tt, with com- 
plex amplitude Ai expi(^z/2, which has been extensively 
used to compare the exact quantum Green's function 
with its semiclassical estimation at the leading order in h. 
In the same way, the Fourier transform of the following 
function: 

ffi(C) = -^^^E^-^c(q)^r,c(qo)r'/'MA-A,(c)) 



-CE-^''^"^^^^' 



; -^(pi) 



(93) 



will also depict peaks at the classical actions Si/2tt, 
whose complex amplitude, given by 



-AiCi exp iq 



(94) 



allows us to extract the numerical value of the h correc- 
tion Gj. 

The energy A being fixed, the 6{X — Ar(C)) function 
selects the values C-r(A) of ( for which A is an eigenvalue, 
transforming Eqs. (p3) and IB3) into 
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ffo(C) 
51 (C) 



4 ^ (r,Cl^k,C> 



(95) 



■E' 



(r,C|i^|r,C> 



C'^/^-J (C - Cr(A)) - CE-^' cos (C^' + M- 



r 



Stepping to the case of the trace of the Green's func- 
tion, the preceding relations (891) and (BG) becomes 



1 



-Im Tr G(A, /i) == V (5 (A - A^(/i)) 

'TT ■*- ^ 



and, see Eq. (|65|): 



ilm TrG'(A, ^i), = — ^ylf (cos ( ^Sf 



^^iC,"' 



1 



(96) 



-Sr + <Pr]}, (97) 



where 5*'' is the action of the periodic orbit and 

To 



A]' = 



|det(m(ro)- 1)1^/2' 



^r 



:m, 



(98) 



L cr = Ci(To) + c/^^(ro) 



so that the classical quantities S\^ ^ A^ , and the h correc- 
tion C\'' can be obtained by taking the Fourier transform 
of the following expressions with respect to the variable 



5r(c) = 



IT Y^ 1 

2^(r,C|^|r,C) 

-IE 



(T,Clf|r,C) 



C^<5(C-C.(A)) 



C'<5(C-Cr(A))-CE-^rcos(C5r 



^D 



(99) 



r 



B. Computing quantum quantities 

Focusing on the A = 2 value, the 2D hydrogen in a 
magnetic field case, one has to find effective h values 
for which 2 is an eigenvalue of the Schrodinger equation 
Tl{h)ip{u,v) = 2ip(u,v), which is conveniently written as 
follows: 



2 + e{u^+v^)--u^v^{u^ + v^) 






'ip{u,v) 

V'(w,w) (100) 



such that a = h^ appears to be a solution of a generalized 
eigenvalue problem {A — aB)ip ~ 0, with 

A = 2 + e{u^ + v^) - -u^v^{u^ + v'^), 

Ifd^ ^x (101) 

The preceding operators A, B, and thus 7Y(/i) are in- 
variant under all transformations belonging to the sym- 
metry group 6*4^, leading to four nondegenerate series 



of energy levels, labeled EEE, EEO, OOE and 000 ac- 
cording to Ref. [^ and a twofold degenerate series EO 
and OE, where E means even and O means odd, the first 
two letters referring to the u — > ~u and v —^ —v sym- 
metries, the third letter to w <-> w. Actually, because 
of the definition of the semiparabolic coordinates {u,v), 
only eigenvectors invariant under the parity symmetry 
ip(—u,-~v) = '0(u,w) correspond to eigenvectors of the 
2D hydrogen in a magnetic field, allowing us, in princi- 
ple, to drop the OE and EO series M, Em. However, from 
the semiclassical point of view, one would have to extend 
all preceding sections to symmetry-projected propagator 
and Green's function |2^, and thus to take into account 
symmetry properties of the classical Green's function, 
which is beyond the scope of this paper. For this reason, 
we also include the OE and EO series in the remainder 
of this paper. 

Finally, eigenvalues and eigenvectors are obtained by 
solving the matrix representation of the generalized 
eigenvalue problem (A — aB)'>jj = in sturmian bases 
(one for each symmetry class) [Q, using the Lanczos al- 
gorithm. Typically, we have computed effective h values 
ranging from to 124, which for scaled energy e = —0.1 
corresponds to roughly 61 000 eigenvalues in total. One 
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must notice that the generaUzed eigenvectors |t, ft), for 
a fixed h value, are actually orthogonal for the scalar 
product defined by operator B = p^/2: 



giving rise to 50,1 (C) ( P5| ) and goa(C) ( pQl ) expressions in 
terms of the computed eigenvectors: 



{r,h\^\r',h)^Srr' 

so that the \t, h) and |r, ft) relations read 

1 



\T,h) 

M) 



|r,ft), 



r, ft|r, ft) 
1 



(102) 



V{r,h\B\T,h) 



(103) 



|r,ft), 



<?o(C) = 

51(C) = 
ffo^(C) - 

gno = 



^^^^^.,c(q)^r,c(qo)C'/'5(C-C.(2)), 



(2^) 



3/2 



^^,,c(q)^-c(qo) c'/' MC - Cr(2)) - cE-^' co«(C^' + '^')' 



5^(r,C|r,C)C'<5(C-C.(2)) 



(104) 



-^^(r,C|r,C)C'<5(C-C.(2))-CE-^*'''^°^(C^r' 



^D 



As explained previously, the Fourier transform of the 
two functions gi and g^ will depicts peaks at classical ac- 
tions and ft corrections are obtained from the amplitude 
of these peaks. However, in the case of signal given by 
c{t) = X] On exp (iujnt), it is now well known that the har- 
monic inversion method is very well suited and is much 
more powerful than the conventional Fourier transform to 
extract unknown frequencies ijJ„ and amplitudes a„ [ pO[ . 
In our case the signals are the two functions gi(C) and 
gf{C), which are of the form J2i^i^i sin ((S"; -I- </>;) be- 
sides contributions from all other types of orbits (ghost, 
continuous family, etc.). 



C. h corrections for G(q, qo,2) 

Orbits having initial and final points at the nucleus 
(i.e., q = qo = 0) are of special interest because they are 
involved in semiclassical estimation of the photoioniza- 
tion cross-section [G5l E9|, which can be directly com- 
pared to experimental results |30, ^. Even if the 



full ft expansion of the cross-section does not reduce to 
G(0, 0, 2) contributions, all closed orbits are well known 
and classified, so that this case remains a nice example 
of ft corrections for G(q, qo, 2). 

The Fourier transforms of both functions 50(C) (up- 
per plot, solid line) and (71(C) (lower plot, solid line), for 



r 



scaled energy e = —0.1, are displayed in Fig. |[ More 
precisely, 50(C) and gi(C) being known only on a finite 
interval [0, Cmax], we have plotted the modulus of their 
windowed Fourier transforms, defined as follows: 



Fois) 



Fiis) 



6 



(Cmax) Jo 



c„ 



dC C(Cmax - C)5o(C)e 



-i2-KsC, 



(Cr --^^ 



max^ Jo 



Cmax 



dCC(Cmax-C)5l(C)e 



-■i27rs(^ 



(105) 



As expected, they depict peaks at the classical actions 
of closed orbits, whose trajectories in (m, v) plane have 
been inserted in the figure, the black circle corresponding 
to the nucleus position. In the figure, the dotted lines 
corresponds to the semiclassical estimations of the same 
functions using the classical properties given by Table |. 
The closed orbits being either half of a periodic orbit or a 
periodic orbit, we label a given close orbit with the four- 
disk code of the corresponding periodic orbit |^, ^ . 

For the leading order in ft (upper plot), as expected, 
the agreement between the quantum results and the semi- 
classical estimation is excellent. For the first order ft 
correction, the agreement is very good, but one can no- 
tice that there is a discrepancy for the amplitude of the 
last two peaks. This is not due to errors or inacurrate 
calculations in the semiclassical estimation, but rather a 
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1.5 

action/2ji 



Table I: Classical properties of closed orbits involved in 
the semiclassical expansion of the quantum Green's function 
G(q, qo, 2) of the 2D hydrogen atom in a magnetic field, for 
the case q = qo = 0. Because each closed orbit corresponds 
either to a half-periodic orbit or a periodic orbit, we have 
labeled them with the four-disk code of the corresponding 
periodic orbit |34^p3[. Their trajectories in the {u,v) plane 
are shown in Fig!^. 5"; is the reduced action (i.e., J p dq/27r), 
Ti is the period, Ai is the leading semiclassical amplitude, vi 
is the Maslov index, Ci is the first order h correction, given 
by the sum Ci(0, 0,Ti) + Cf""-^(0, 0, T;), see Eq. (|l|). 



Code 



3i 



Ti 



A, 



13 



1243 



1.094 570 5 

1.564 998 2 

121343 1.791060 7 

12124343 1.933 5221 
1212134343 2.031948 2 



2.425 093 3 
3.6001374 
4.286 257 7 
4.796 775 8 
5.214 323 3 



I'l 



0.295 342 6 1 

0.152 365 2 

0.109 503 9 3 

0.093 368 7 4 

0.086142 5 



Figure 3: Modulus of the windowed Fourier transforms Fo 
(solid line, upper plot) and Fi (solid line, lower plot), see 
Eq. (105), of the quantum functions go (lea ding order in h) 
and gi (first order h correction), see Eq. ([104), associated 
with the quantum Green function G(q, qo,2) in the case of 
the 2D hydrogen atom in a magnetic field and for q = qo = 
(see Sec. IVJ for all details). As expected from semiclassical 
formula (eBD, peaks are appearing at action (i.e. Jpdq/27r) 
corresponding to classical orbits having initial and final posi- 
tions at the nucleus. For the first five ones, the trajectory in 
the {u, v) plane are also plotted, the nucleus being depicted 
by the black circle. The agreement with the semiclassical es- 
timations of these functions (dotted lines) is excellent, even if 
discrepancies in the amplitude of last two peaks in the lower 
plot can be observed. These are actually a manifestation of 
limitation of the Fourier transform and not inacurrate calcu- 
lations of the h corrections, as it is emphasized by the quan- 
titative comparison (using harmonic inversion) displayed by 
Table 0. 



manifestation of the limitations of the Fourier transform. 
To emphasize this point, we have used the harmonic in- 
version to extract, for each of these orbits, the h correc- 
tion coefficients C^^ , from the quantum function (?i(C)- 
The results are compared to the classical calculation Ci 
in Table O. The agreement is excellent, the relative er- 
ror on the amplitude being lower than 10^^. As usual, 
the phase extracted using harmonic inversion, being the 
most sensitive quantity, the agreement on the sign of the 
Ci, rather nice for the first four orbits, decreases rapidly. 
Finally, one must mention that this good agreement be- 
tween quantum and semiclassical calculations has also 
been found when considering quantum Green's functions 
G(q, qo,2) with other initial or final points. 



Code Ci(o,o,ro Ci^--^(o,o,ro 



C, 



13 



-0.202 769 9 
-0.119 409 3 
-0.148 282 2 

-0.172 990 6 

1212134343 -0.192 904 3 



1243 

121343 

12124343 



0.016 539 4 
0.019 7412 
0.041175 5 
0.071748 
0.1174645 



-0.186 230 5 
-0.099 6681 
-0.107106 7 
-0.101242 7 
-0.075 439 8 



Table II: Numerical comparison between the theoretical h cor- 
rections Ci for the quantum Green's function G'(q, qo, 2) of the 
2D hydrogen atom in a magnetic field, for the case q = qo = 
and the numerical coefficients Cf^^ extracted from exact quan- 



tum function (?i(C) (Eq. (104)) using harmonic inversion (tak- 



ing into account multiplicity) . The agreement is excellent for 
the amplitudes and rather nice on the phases, thus empha- 
sizing the validity of the present theory. That the agreement 
becomes less good for the last orbit only shows the limitations 
of the harmonic inversion method, which usually appear on 
the phase. 



Code 


Ci 


icn 


Rel. error 


argQ^^ 


13 


-0.186 230 5 


0.1864 


f»8 X 10"-* 


1.002x7r 


1243 


-0.099 6681 


0.0995 


fa 2 X 10"^ 


l.OlXTT 


121343 


-0.107106 7 


0.1072 


f« 9 X 10"^ 


1.02x7r 


12124343 


-0.101242 7 


0.1016 


f»4 X 10~^ 


1.04x7r 


1212134343 


-0.075 439 8 


0.0761 


Ri9 X 10"^ 


1.14x7r 



plot, solid line), defined, as previously, as follows: 



Fl'is) 



6 



(C„ 



(C„ 



dC C(Cmax - 


-C)5o''(C)e-*'"^^, 


dC C(Cmax - 


-osrxoe-*'"^^. 




(106) 



D. 



h corrections for Tr G(q, q, 2) 



Still working at scaled energy e = —0.1, Fig. depicts 
the modulus of the windowed Fourier transforms of gQ^{C) 
and g\'^{C), Fq" (upper plot, solid line), and F*'' (lower 



The trajectories in the (w, v) plane associated with the 
peaks are also plotted in the figure. The classical prop- 
erties of the corresponding periodic orbits are displayed 
by Tabic III, Again the agreement is excellent between 



the quantum results (solid lines) and the semiclassical 
estimation (dotted lines). The quantitative comparison 
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3.0 

action/27t 



Figure 4: Modulus of the windowed Fourier transforms Fg^ 
(solid line, upper plot) and Fi"^ (solid line, lower plot), see 
Eq. (lOf)), of the quantum functions Qq^ (leading; order in h) 
and gi (first order h correction), see Eq. ( |l04[ ), associated 
with the trace of the quantum Green's function Tr G(q, q, 2) 
in the case of the 2D hydrogen atom in a magnetic field (see 
Sec. [Vll fo r all details). As expected from semiclassical for- 
mula (pq), peaks are appearing at action (i.e. ^pdq/27r) 
corresponding to classical periodic orbits, whose trajectories 
in the (u, v) plane are plotted ( the nucleus being depicted by 
the black circle). The agreement with the semiclassical esti- 
mation (dotted lines) is excellent, as it is emphasized by the 
quantitative comparison (using harmonic inversion) displayed 
by Table fvl. 



and the values C^^^ 
extracted from the quantum function 5*'^(C) 



between the classical coefficients C*"" 



Table [V 



IS given m 
The agreement is excellent for the amplitude 
of the coefficients and is rather good for their phases, 
which emphasized the validity of the semiclassical for- 
mula developed in the preceding sections, especially the 
additional term arising from the Jacobian describing the 
change from the Cartesian to local (along the periodic 
orbit) coordinates (see Eq. (p8|)) and which contributes 
to a large part of the ft correction for the present orbits. 



VII. CONCLUSION 

In summary, we have explained in this paper how to 
effectively compute ft corrections in the semiclassical ex- 
pansions of the propagator iir(q, qo,T), its trace K{T), 
the quantum Green's function G'(q, qo,i?) and its trace 
G{E) for chaotic systems with smooth potential. The 
method is based on the classical Green's functions asso- 
ciated to the relevant trajectories, that is either going 
from q to qo in the propagator case or periodic orbits for 
K{T), together with adapted boundary conditions. We 
have shown how all quantities can be obtained by inte- 
grating, using the standard Runge-Kutta method, sets 



Table III: Classical properties of periodic orbits involved 
in the semiclassical expansion of the trace of the quantum 
Green's function Tr G(q, q, 2) of the 2D hydrogen atom in 
a magnetic field. Their trajectories in the (it, v) plane are 
shown in Fig. H. Sj^ is the reduced action (i.e,. ^pdq/27r), 
Ti is the period, AY is the leading semiclassical amplitude, fj,i 
is the Maslov index, C;'" is the first order h correction, given 
by the sum Ci{Ti) + CT^'^iTi), see Eq. (E|). 



Code 



sr 



ytr 



AT 



Mi 



1234 


2.709 8513 


6.204155 6 


0.827 8814 


4 


1243 


3.129 996 4 


7.200 274 7 


0.616 496 8 


4 


12434 


3.2271681 


7.541640 6 


0.548 4791 


5 


123434 


3.272 2381 


7.748 406 8 


0.555 880 6 
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Code 


Ci{Ti) 


Cf^^(TO 


cr 




1234 


-0.622577 


0.026912 


-0.595665 




1243 


0.166821 


0.051665 


0.218486 




12434 


-0.203536 


0.058541 


-0.144995 




123434 


-1.41705 


0.07241 


-1.34464 





Table IV: Numerical comparison between the theoretical h 
corrections C;'' for the trace of the quantum Green's function 
Tr G{q, q, 2) of the 2D hydrogen atom in a magnetic field and 
the numerical coeffici ents C/^^ extracted from exact quantum 



function gi{C) (Eq. (104)) using harmonic inversion (taking 
into account multiplicity) . The agreement is excellent for the 
amplitudes and rather nice on the phases, thus emphasizing 
the validity of the present theory, especially the additional 
term due to the transformation from the Cartesian coordi- 
nates to the local frame along the periodic orbit (see Eq. (kq)). 



Code 


cr 


\cr\ 


Rel. error 


argC,^^ 


1234 


-0.595 665 


0.5958 


« 2 X 10-"* 


1.005x7r 


1243 


0.218 486 


0.2178 


^ 3 X 10-=* 


0.04x7r 


12434 


-0.144 995 


0.147 


^ 1 X 10"^ 


0.93x7r 


123434 


-1.344 64 


1.347 


^ 2 X 10"^ 


O.OSXTT 



of differential equations. We have also shown that in the 
derivation of the semiclassical expansion for K{T) (and 
thus G{E)), starting from the Feynman path integral, 
one must take into account additional terms, which af- 
fect only ft correction coefficients. This is emphasized by 
the excellent agreement observed when comparing, in the 
case of the 2D hydrogen atom in a magnetic field, our the- 
oretical results with the numerical coefficients extracted 
from exact quantum data, using the harmonic inversion. 
Obviously, there are still many points to be developed. 
Besides the few cases, such as self-retracing orbits or con- 
tinuous families of orbits, needing specific extensions, it 
would be very interesting to understand how to include 
continuous and discrete symmetries. Also, going into the 
extended phase space (q, t, p, —E) ||2j , it would be possi- 
ble to get a better understanding of similarities observed 
between the differential sets leading, on one side to the ft 
corrections for the propagator and its trace and, on the 
other side to the additional terms arising in the ft correc- 
tions for the quantum Green's function and its trace. 
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In the case a\\ = 0, we have thus found another eigenvec- 
tor for the eigenvalue 1, which means that a small dis- 
placement of initial conditions in the ej^ direction leads 
to another periodic motion with the same period T, and 
thus that the periodic orbit is actually embedded in a 
continuous family. Indeed, using notations from Sec. |y|, 
one can show that 



XW(0) = -Me^ 
"II 



(A5) 



so that we have 



Appendix A: FEW PROPERTIES OF M{T) 

In this appendix, we consider an isolated unstable peri- 
odic orbit of period T. We shall use the notations e|| and 
e_L for the units vectors which are, respectively, parallel 
to the flow and perpendicular to the energy shell at the 
initial point. From Hamilton's equations, we have that 
M(T) ■ e|| = e||, i.e., e|| is an eigenvector of the matrix 
M{T) for the eigenvalue 1. The symplectic equation ful- 
filled by M{T), namely, Af (T)^ • S • M{T) = E, implies 
that, if e^ and Gj are two eigenvectors for the eigenvalues 
\i and Aj, we have the following properties: 



M{tY ■ (SeO = -^ (Se,) 
(A,A, - l)e7Se, = 0, 



(Al) 



showing thus that 1/Ai is an eigenvalue of Af (T)^ and, 
from that, of Af(T). In addition, M{T) being a real 
matrix, \j and 1/Aj are also eigenvalues of M{T), so 
that the nontrivial eigenvalues (i.e., y^ 1) either fall in 
the (A, 1/A) pair or in quadruplet (A, 1/A, A, 1/A). 

In the case of e^ — e||, the two preceding equa- 



tions ( |Al| ) imply that ej_ is an eigenvector of M{T) 
(but not necessarily of M{T)) for the eigenvalue 1 and 
that for every Aj 7^ 1, Gj is orthogonal to ej^. In the basis 
(e||,e^,ei,- • • ,e2/-2), M{T) entries then read 



M{T) 



1 a|l 

10 

ai Al 

^2 A2 







a2/-2 







A 



2/-2 



(A2) 



where we have supposed that all eigenvalues are simple. 
For degenerated eigenvalues, M{T) would be block diag- 
onal. For a generic periodic orbit, an and ai are nonvan- 
ishing emphasizing thus that e_L is not an eigenvector of 
M{T). Introducing the vector e_L defined as follows: 



e± 



2/-2 

E 



Ps^i 



with 



/3. 



1- A, 



one immediately gets that 

M{T)e± = e± 



aiieii 



(A3) 



(A4) 



an 



\X\fdET. 



(A6) 



In Sec. M, one needs to compute derivatives with re- 
spect to the period T of det (rn[T) — 1), whose expression 
in terms of the nontrivial eigenvalues of the monodromy 
matrix reads 



2/-2 

det (m(r) - 1) - n (^J- - !)• 



(A7) 



Introducing Vn and 'P_\_ the projectors on the directions 
e|| and ej^, more precisely 



V\\ 



and 



Ti 



(A8) 



'II =^11 -en 

ones defines the matrix N{T) as follows: 

N(T) = M{T) -{l-V\\-Vi_). (A9) 

In the basis (e||, ej^, ei, • • • , e2/-2)j using orthogonal- 
ity between e^ and e^, entries of N{T) read 



N{T) 



1 a\\ 71 72 

10 

ai Al - 1 

a2 A2 - 1 



a2/-2 



72/-2 







A2/-2 



(AlO) 



where 7j = eJ • e^, which actually could be related to 
the Uj, but this is not necessary in our case. This shows 
that the determinant of N{T) is exactly JljiT ('^i ~ !)• 
The main ad vantage of the matrix N{T) is that its ex- 
pression (A9) does not involve the eigenvectors or the 
eigenvalues of M{T), so that its determinant can be di- 
rectly computed, without the diagonalization stage re- 
quired when getting det {m(T) — 1) through the eigen- 
values Xj. Furthermore, derivatives of IndetA^(T) with 
respect to the period T are also straightforward to ob- 
tain, knowing derivatives of M{T) and of X(T), whereas 
derivatives of Xj would require the knowledge of those of 
the eigenvectors e^. 
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